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Investigation of the Temperature Dependence of Axial Dechanneling 
by Point Defects 


Calculations of the depth dependence of the Rutherford back- 
scattering minimum yield were performed for crystals with 
non-correlated displaced lattice atoms (point defects). The 
main result is that the slope of the temperature dependence 
changes considerably for different displacement distances ae 
In order to explain this effect the contributions of the 
different dechanneling processes to the minimum yield are 
discussed in detail. It is shown that the combined action of 
these processes is quite different for different values of ra 
which is responsible for the strong influence of the value of 
r, on the temperature dependence of the minimum yield. This 
effect can be used for the analysis of point defects with 
preferred positions in the lattice cell. For example defects 
in N implanted GaAs are analysed. 
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4. Introduction 


For about 20 years channeling has been used fou detect stu- 
dies in crystals [1-9]. The relative Rutherford backscat- 
tering yield for ions incident parallel to a crystal axis 
CRBS minimum yield Kan is increased if the crystal contains 
defects (displaced lattice atoms). The enhancement of the 
minimum yield is caused by direct backscattering and dechan- 
neling of the ions due to the interaction with the displaced 
lattice atans. The dechanneling mechanisms are different for 
different kinds of defects which are characterized by diffe- 
rent correlations between the displacements of the lattice 
atoms. This paper deals with uncorrelated displaced lattice 
atans (point defects, small defect clusters). With respect to 
dechanneling point defects (and small clusters) are characte- 
rized by the relative number nyd of lattice atoms which are 
displaced from their lattice sites and by their displacement 
distance r, perpendicular to the channeling direction (as- 
sumption of cylindrical symmetry is mostly justified). Only 
in the case where Pr. can be assumed to be randomly distribu- 
ted within the area belonging to one atomic string (heavily 
damaged crystals) the defect density yal) can be calculated 
directly from the minimum yield X ant®) measured as a func- 
tion of the depth z [1,2.3,7]. Otherwise, additional informa- 
tion is necessary for the determination of both the defect 
density Nyglz) and the distribution of the displacement 
distances. rae Diffusion model calculations of Matsunami et 
al. [10,11] proved the temperature and energy dependence of 
dechanneling to depend on the value of r,- This offered the 
possibility of a quantitative defect analysis by performing 
RBS measurements at different temperatures and different 
energies. More detailed investigations using a discontinuous 
model [7] showed that- only the temperature dependence of 
dechanneling provides a useful information about the distri- 
bution of r,. Therefore, the temperature dependence of axial 
dechanneling in crystals with non-correlated displaced lat- 
tice atoms is of special interest in this paper. The investi- 
gations are performed assuming all displaced lattice atoms to 
have the same value of the displacement distance ra Gar irs 
distributed). This corresponds to point defects with. prefer— 
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red positions in the lattice cell which appear mainly in 
weakly damaged crystals. The results are discussed and ap- 
plied to the analysis of defects in Si and Gads. 


2. Dechanneling model 


The model used is described in ref.{7] and in the following 
it is referred to as the discontinuous model. It is based on 
the general description of dechanneling given by Lindhard 
[12]. Instead of a continuous distribution g(Ey,z) of the 
transverse energy ŒE, of the ions at depth z the relative 
number g,(z) of ions with transverse energy between Ey,-4 and 
Bay (Ey, = EL, 1/ G ts io ax? see'ref.[7]) is the 
physical quantity of interest. It is determined by a system 
of differential equations (master equation) 


dalal a Saeed, a) 


where the dechanneling matrix Qij is given by the probability 
per depth Pag for the transition of an ion from group j to 
group i 


The Rutherford backscattering minimum yield X is obtained, 
by 
X minlz) =È Te (2), (3) 


where T, is the relative probability for hitting a regular 
lattice atom or a displaced atom [7]. The dechanneling matrix 
is the sum of all contributions from statistically indepen- 
dent processes (electronic scattering, nuclear scattering, 
scattering on point defects) 


Qs = 0455°? + Crm gC2d]Q, + ni g(ar0, P(r.) (4) 


The calculation of ig 2, PAN and agg Pe Gr) is descri- 
bed in ref.[7]. Note that the dechanneling matrix contribu- 
tion of the point defects depends on the displacement distan- 
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3. Temperature dependence of dechanneling 


The temperature dependence of dechanneling is quite different 
for different values of r,: This is demonstrated in fig.1a 
for the example of 1 MeV H* channeling in <111>Si with 10% of 
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Fig.1. Temperature dependence of Xin (fig.1a) and AX oon 
(fig.1b) at depth z = 2, um for 1 MeV H? incident on perfect 
<411> Si and <111> Si with different point defects. 


the lattice atans being displaced by rą = 0.025 nm and 


ie 0.055 nm, respectively. Of course, in all cases the 
minimum yield Xan increases with respect to the minimum 
yield for the perfect .crystal Capa = 0). However, the tem- 
perature dependence of Xain is weaker or stronger than that 
of the perfect crystal Tomir = 07025 nm or - rodz 0.055 nm, 
respectively. This can be seen more pronounced in fig.1b 


which shows the difference of the minimum yields 


ax x 


min = “min FR ae pe rv The slope of AXE depends consi- 
derably on the value of r,, and it changes even its sign. An 
explanation of this effect is given in the following. 
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The dechanneling contributions due to the scattering of the 
ions at the electrons (e), at the thermally vibrating lattice 
atas (nuclear scattering, (n)), and at the displaced lattice 
atoms (point defects, (d)) are demonstrated in fig.2. 


transverse energy distribution gE, arbitrary units) 


transverse energy E, lev! 


Fig.2.) Transverse energy distributions at depth z = 0 and 
2 = jum for 1 MeV H} incident on perfect <111> Si (fig.2a) 
and on <111> Si with 10% lattice atoms being displaced by 
ra = 0.025 nm (fig. 2b) and by r, = 0.055 nm (fig.2c). The 
temperature is 160K in all cases. For the demonstration of 


the different dechanneling processes by the arrows see the 
text. q 


In the transverse energy concept dechanneling means the 
change of Ey by the processes mentioned above. For a qualita- 
tive discussion the minimum yield can be estimated by 
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oo 

Xan) Ae fai eC 84,2), (5) 

Erc 

where Ey, = 3 Eye is the critical transverse energy CY, is 
the critical angle introduced by Lindhard}. 
The initial transverse energy distribution (z = 0) is heavily 
peaked at small transverse energies. Due to dechanneling the 
transverse energy distribution changes with increasing depth. 
However, small Ey, are still dominating. The situation for 
2 = 4 um is shown in fig.2. It can be seen that in all three 
cases the enhancement of the minimum yield with increasing 
depth is caused by the reduction of the peak of Bl Ey ) at 
small E}. Therefore, with respect to dechanneling the main 
effect is the transverse energy transition from FE, << Ete to 
Ey > Ey. This is provided by the three processes Ge) an), 
and (d) mentioned above. However, the relative contributions 
of these processes are different for different transverse 
energy regions and for different values of ro. This is demon- 
strated in fig.2 by arrows of different positions, lengths, 
and thicknesses. i 
(i) Perfect crystal: For small E, the transverse energy is 
changed mainly by electronic scattering. The contribution of 
nuclear scattering is less by about two orders of magnitude 
(ions with small E do not reach the thermally vibrating 
lattice atoms). The changes of By due to electronic scatte- 
ring are small. Therefore, the electronic processes are 
effectiv only in a snall E, region up to medium transverse 
energies (see arrow (e) in fig.2a). With increasing Ey the 
nuclear processes doaninate more and more. Their effective 
Ey region (arrow (n,T) in fig.2a) increases to lower Ey for 
increasing temperature. In summary, the enhancement of the 
minimum yield Xin is the result of the combined action of 
the electronic and nuclear scattering processes. The tempera- 
ture dependence of Ye is introduced by the nuclear proces- 
ses. 
(ii) Crystal with slightly displaced lattice atoms (e.g. 
rą = 0.025 nm): In comparison with the perfect crystal (i) 
there are additional scattering processes at the displaced 
lattice atoms. For small displacement distances GE 2 thermal 
vibration amplitude) the effective Ey region for ion scatte- 
ring at the displaced lattice atoms is similar to that for 
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nuclear scattering. With increasing r, this region extends to 
lower E, (arrow (d) in fig.2b). The influence of the addi- 
tional thermal vibration of the displaced lattice atoms on 
the dechanneling was shown to be negligible for 
Ta = 0.025 nm. Therefore, the scattering at these defects is 
approximately independent of the temperature. For small By 
the transverse energy is changed only by electronic processes 
as in the case of the perfect crystal. The further enhance- 
ment of the transverse energy up to Ey > EL, is provided by 
the nuclear scattering (n) and by the scattering at the de- 
fects (d). As it is expressed by the thicknesses of the 
arrows the resulting minimum yield in fig.2b is larger than 
that of the perfect crystal (fig.2a), but its contribution 
(n) from the scattering at the thermally vibrating lattice 
atans is less than in the case of ‘the perfect crystal. There- 
fore, the temperature dependence of Lr for crystals with 
slightly displaced lattice atoms is weaker than that for the 
perfect crystal as shown in fig.1 (negative temperature 
dependence). 

(iii) Crystals with heavily displaced lattice atoms (e.g. 
Taar OROS nm): Defects with large displacement distances 
r, % 0.05 - 0.10 nm interact preferentially with ions of 
small transverse energy. Therefore, their effective Ey region 
is similar to that of electronic scattering (fig.2c). This 
causes much more ions with small Ey to increase their trans- 
verse energy than in the case of the perfect crystal. The 
further enhancement of the transverse energy up to Ey, > EL. 
is obtained only by nuclear processes. Because of the large 
transverse energy transfer for small Ey the nuclear processes 
are much more effective than they are in the perfect crystal. 
Therefore, in this case Xin is much larger and the tempera- 
ture dependence of Mee A is stronger than for the perfect 
crystal according to the results. depicted in fig.1. 


4. Application to defect analysis 


As shown in chapter 3 the temperature dependence of 
AXnin > Tinie E Xmin, perf 
values of r„. Therefore, both Da qb?) and r, can be determined 


is quite different for different 


a 
from AX va nl? T) measured caus different temperatures. The 
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procedure is described in the following. 

In the first step the difference of the minimum yields 
BM oz slp 2 is measured for a fixed temperature Ths Assuming 
a set of different values of T, the corresponding defect 
densities ng bz) are calculated from Sa Cao The result 
is a-set of (z,, nala). In the second step AX inte?) is 
calculated for different temperatures using the set of (ra? 
npg lZ). The calculated temperature dependence of aX vin is 
quite different for different Cres n a62))- The comparison 
with Po BME.) measured for T T, provides the correct r 
and nalz). In order to get sufficient accuracy measurements 
should be performed in a temperature interval of 100K at 
least. 

In fig.3 and 4 this method is demonstrated for the example of 
He implanted silicon using Tp, = 160K. The calculated tempera- 
ture dependences of AX iin at- depth z = epum for different 
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Fig.3. Temperature dependence of the calculated BET at 
depth z = 2,um for different values of Da (Galen SEU! ebe\sl 
wee at depth z = 2,;umsas a function of T (fig.3b) for 
1 MeV H° incident on He implanted <111> silicon. The calcula- 


ted curve aX. kr.) contains the statistical error introdu- 
ced by simulation. For the experimental data see fig.4. 
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Fig.4. Depth dependence of the measured and calculated AX en 
at temperatures of 40K and 280K for 1 MeV H* incident on 
<111> silicon implanted with 2 MeV He ions (3.6#10' Com, 
40K) and annealed at 290K for 10 min (fig.4a) and the depth 
dependence of the point defect density Noa calculated for 


r, = 0.055 nm (Eig Ub) 


values of r, (and the corresponding n qo), not shown) are 
depicted in fig.3a. The comparison ith the experimental data 
of Howe et al.[13] provides r, = 0.55 nm, For the estimation 
of the error the second difference ax 4 = AX C2 jum, 
280K) - AX, (Zum, 40K) as a measure of the slope of the 
temperature dependence of AX win is amore useful quantity 
(see fig.3b). It can be seen that Noo a depends heavily on 
Tą for r; 5 0.07 nm. For r, > 0.07 nm AA in changes weakly 
and is approximately equal to the value obtained for ra 
randomly distributed within the area belonging to one atomic 
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string. Therefore, the method considered is very sensitive 
for point defects with small and medium displacement 
distances (rą < 0.07 nm for <111>Si). Fig.4 ‘shows the 
functions nig(z), AX 62 40K), and AX. 6 280K) calculated 


for r, = 0.055 nm and the corresponding experimental data. ` 


The Aad agreement for the two temperatures is obtained only 
for this value of r,» 

As a second example point defects in <100>GaAs implanted with 
280 keV N at room temperature are investigated for different 
nitrogen doses. The results obtained for 7440! 3om7? N implan- 


tation are depicted in fig.5a and 6 (Tp = 295K). 
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Fig.5. Temperature dependence of (et calculated for diffe- 
rent values of r, measured with 1.4 MeV He” incident on 
<100>GaAs implanted at room temperature with ™ 401 3Nom72 
(fig.5a) and with 10! Nom? (GPa aoao 


They proof the existence of a negative temperature dependence 
of oe The conparison of the measured and calculated 
temperature dependence of Xi provides ae 0.022 nm. The 
corresponding calculated npa EA and axis nt? ,I) are given in 


fig.6. For the higher dose of 10! 8Nom~ az the results (figs.5b 
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Fig.6. Depth dependence of the measured and calculated BE in 
at temperatures of 125K and 295K for 1.4 MeV Het incident on 
<100>GaAs implanted at room temperature with 7101 Sim 
(fig.6a) and the depth dependence of npg calculated for 
r, = 0.022 nm (fig.6b). 


and 7, Tp = 295K) show a positive temperature dependence of 
Xun providing r 0.065 nm. The corresponding calculated 
ng 62) and AKinl zD are depicted in fig.7. These results 
indicate that the N implantations with different doses cause 
different point defect complexes which are characterized by 
different displacement distancés r.. Possible defect comple- 


a 
xes are discussed in ref.(14]. 
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Fig.7?. Depth dependence of the measured and calculated aX in 
at temperatures of 125K and 295K for 1.4 MeV He* incident on 
<100>GaAs implanted at room temperature with 101 oN om? 
(fig.7a) and the depth dependence of Noa calculated for 
oe ge 0.065 nm (fig.7b). 
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